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In this Letter we show that the evidence for neutrino oscillations from the Super- 
Kamiokande atmospheric neutrino data fully determines the 3x3 neutrino-oscillations 
mixing matrix and predicts an energy independent solar neutrino deficit at the level of 
45%. This corresponds to a ratio of measured to predicted neutrino flux of R^ olar = 0.55, 
in good agreement with the experimental results. We achieve this result within the frame- 
work of a minimal, three-generations neutrino mixing, with mass squared differences of 
AM 2 ~ 0.45 eV 2 and Am 2 = 0(1O -3 ) eV 2 . The mixing matrix derived here is char- 
acterized by the mixing angles 8 = 35.1°, j3 = 5.5°, and ip = 23.3°, and a vanishing 
CP-violating phase, 5 = 0. 
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1. Introduction 

The phenomenon of neutrino oscillations was first postulated by PontecorvoEl in 
f957, and has been ever since a prime candidataJn explaining the long-standing 
puzzle provided by the solar neutrino deficit .BooEm The anomalies reported in the 
atmospheric neutrino production0€na and the positive sig nals obtained by LSND in 
both the Va — > D e and the — > v e appearance channelsEJJO have contributed to fur- 
ther strengthen the evidence towards neutrino oscillations. At the same time, these 
observations raise the question whether they can all be explained in the minimal, 
three-generation neutrino mixing formalism. Typically, all experiments have been 
analyzed in terms of the simpler, two-neutrino oscillation hypothesis, which in turn 
have given rise to three different mass squared differences, Sm 2 . The solar neutrino 
deficit implies a Sm 2 of 0(1O~ 10 ) eV 2 or (D(10~ 5 ) eV 2 , corresponding to vacuum or 
MSW oscillations, respectively. The atmospheric neutrino deficit indicates a dm 2 
of 0(l(r 3 ) eV 2 . The LSND evidence points towards a Sm 2 of 0(1) eV 2 . Since only 
two mass squared differences are independent, one possible solution to this puzzle is 
the introduction of a fourth neutrino flavor, which must be sterile if its mass is low. 
Another solution is to ignore one of the indications listed above, typically LSND, 
as it is the only evidence that has not been independently confirmed by another ex- 
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periment. (However, one should keep in mind that the LSND evidence comes from 
two different channels, with different neutrino fluxes, signatures, backgrounds and 
systematics.) On the other hand, if one assumes that the LSND excess is indeed 
due to neutrino oscillations, the two mass squared differences implied by it and 
by the Super-Kamiokande atmospheric data would necessarily lead to an energy 
independent oscillation probability for the solar neutrinos. _____ 

The latter approach has been followed by several authorsJ__l_J__ Using three 
data points, (1) a flat solar neutrino ratio of Rf olar = 0.5, (2) the atmospheric 

deficit for up-going muons (interpreted as —> v r oscillations), and (3) the 
oscillation probability reported by LSND, they fully determine the neutrino-mixing 
matrix. This, together with the LSND and Super-Kamiokande mass squared differ- 
ences, is shown to be consistent with all current neutrino oscillations experiments. 
Alternatively, a more elaborate fit to all datajil yields similar results. What we 
show in this Letter is that the Super-Kamiokande atmospheric data alone contains 
enough information, through its explicit L/E dependence, to fully determine the 
neutrino-oscillations mixing matrix, and is able to predict the solar neutrino deficit 
to a level which is consistent will the current measurements. 



2. Formalism 

Throughout this paper we shall assume that the three neutrino weak eigenstates 
(y e , v^Vt) are linear superpositions of three mass eigenstates {v\, 1^2,^3), 

3 

\v t >=Y J Ut 3 \v j > {l = e,n,r), (1) 

with masses m\ < mi < 771.3. The unitary matrix U, parameterized in terms of 
three mixing angles (0,/3,i/ ) ) and a CP-violating phase (S), reads 

CgC/3 SgCp Sff 

U = ( -CgSf3Sipe lS - SgCip CgC^ - e lS SgSpS^ CpS^e 15 | (2) 

-cgspCjp + ses^e~ lS -sgspc^ - ces^e~ lS cpc^ 

in the Maiani representation,-! where eg — cos 8, sm_^— sm9, etc. Consequently, 
assuming that the mass eigenstates are relativistic_3 and stable, the oscillation 
probability from a state \vi > to a state \i/i> > is given by 



i=l j=i-j-X 



2//, ( ir>;f7,r ( ' ( r,,)M.. ( 2.54 Am?. | 



(3) 



In Eq.(g) above E is the neutrino energy (expressed in MeV), L is the distance 
between the generation point and the detection point (expressed in m), and Am?- = 
mj — m 2 A (expressed in eV 2 ). The neutrino mass spectrum is characterized by two 
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mass squared differences, which we shall loosely refer to as "mass scales" henceforth, 
Am 2 and AM 2 : 



Am 2 .! 



TO5 



Am 2 

AM 2 ee Am^ =m%-m\ 



2 - "n> 

2 _ A 2 _ 2 2 



(4) 
(5) 



The Am 2 mass scale is taken to be of 0(1O~ 3 ) eV 2 , as suggested by the Supcr- 
Kamiokande atmospheric data, while the AM 2 mass scale is taken to be of 0(1) eV 2 , 
as suggested by LSND. Furthermore, we shall only consider the case in which there 
is no CP-violation in the neutrino sector, and thus we shall set the phase 5 = 0. 

3. The Super-Kamiokande Atmospheric Neutrino Data 

Let us assume that at the top of the atmosphere, at t = 0, the number of muon and 
electron neutrinos produced are iV M and N e , respectively. Although both neutrinos 
and antineutrinos are produced in both flavors, we shall use the terms "muon neutri- 
nos" and "electron neutrinos" loosely, to include both neutrinos and antineutrinos. 
The ratio of muon to electron neutrinos, r, 



N(v e + v e ) N e 



(6) 



is, in general, a monotonically increasing function of energy. However, for the 
relevant range of L/E values in Super-Kamiokande, r may be assumed constant 
(r = 2), as we shall do throughout this Letter. Within a detector at a distance 
L ~ t from the production point, the number of muon and electron neutrinos, 
and N' el respectively, are given by: 



N' 



N e P eu , 



N' e = N e P ee + N^ e , 
where P w is the survival probability: 

L 



up, 



4C/ 2 1 C/ 2 2 sin 2 



4 U%U% sin 2 



1.27 Am 2 



E 



4[/ 2 1 C7 2 3 sin 2 



1.27 (AM 2 - Am 2 ) 



P efl is the neutrino oscillation probability P{v e 



E 



efj. 



^U el U e2 U^U^ sin 2 ^1.27 Am 2 



^UeiUesU^U^ sin 2 1.27AM 2 



4J7 e2 C/ e3 C/ A1 2^3 sin 2 



(7) 
(8) 



1.27 AM 2 



(9) 



1.27 (AM 2 



Am 2 ) | 



(10) 
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and P ee is the v e survival probability: 

Pee = 1 - 4C/ e 2 iC/ 2 2 sin 2 N.27Am 2 I) - 4£/ 2 1 [/ 2 3 sin 2 M.27 AM 2 ^ 



4(7 e 2 2 C/ 2 3 sin 2 



1.27 (Ail/ 2 - Am 2 ) — 
E 



(11) 



Since we are only considering the scenario with a vanishing CP-violating phase 5 = 
0, P^e = P e /i- The ratios of measured to predicted (without neutrino oscillations) 
/^-like and e-like events at Super-Kamiokande, i? M and -R e , respectively, are 

K 

R^ = j^ = P^ + 0.5F eAI , (12) 
AT' 

R e = 1 f = Pee + 2.0P^, (13) 

where we have explicitly set r — 2. It is precisely these two ratios, and in particular 
their L/E dependence - as measured by the Super-Kamiokande group - that wc 
shall exploit in order to fully determine the neutrino oscillations mixing matrix. 

For the mass scales considered here, something remarkable happens in the 3rd 
and in the 8th bin of the Super-Kamiokande atmospheric neutrino data - see Fig. 4 
in Rcf.[9]. In the 3rd bin, corresponding to 10 15 m/McV < L/E < 10 2 m/MeV, 
the oscillations involving the Am 2 mass scale are approximately zero, and thus 
negligible, whereas the oscillations involving the AM 2 mass scale average out to 
1/2. One can certainly argue about how many oscillations it would take in a given 
L/E interval in order to safely assume that the oscillations do indeed average out. 
For instance, for AM 2 = 0.18 eV 2 one has approximately five oscillations in this 
interval, whereas AM 2 = 0.36 eV 2 yields obviously twice as many - approximately 
ten oscillations. In our view, even these relatively low numbers of oscillations do 
allow to effectively average out the sine squared variation. Therefore, we shall 
henceforth impose a weak lower bound of about 0.27 eV 2 for the AM 2 mass scale, 
as the average of the two values mentioned above. As we shall see later on, the 
value of AM 2 favored by this analysis lies well above this limit. In the 8th bin, 
corresponding to 10 4 m/MeV < L/E < 10 4 5 m/MeV, all oscillations average out 
to 1/2. (For example, for Am 2 = 10 -3 eV 2 one has approximately nine oscillations 
in this L/E interval, and thus our assumption is well justified.) 

Following these arguments, the ratios of measured to expected /i-like and e-like 
events, as defined by Eqs.(p^)~(p^), are independent of the underlying mass scales 
in these two bins. For the 3rd bin, these ratios yield a rather simple form: 

= 1 - C 2 S 2 (2-4- 2c 2 s 2 ), (14) 
R& = l - 2c 2 s 2 (1 - 2 S 2 ). (15) 

They depend only on the 1-3 and 2-3 mixing angles, (3 and ip, respectively. One could 
argue that using the values measured by Super-Kamiokande, R^ = 0.908 ± 0.068 
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and Re = 1.225 ± 0.108, directly determine these two angles. However, at this 
point we prefer to only restrict the possible values of f3 and ip from the ratio of 
ratios, /R e 3) = 0.741 ± 0.088. This is a much more robust quantity as the 
large systematic errors in the individual flux calculations largely cancel out. The 
allowed values for the mixing angles (3 and tp are shown in Fig. ffl, at the ±lcr 

(8) 

level. The expressions for the u-like and e-like event ratios in the 8th bin, Rl ' and 
R e , respectively, are rather long and cumbersome, and therefore we choose not to 
display them explicitly. With four data points (the ratios of ratios and the /i-like 
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Fig. 1. Allowed regions for the mixing angles (3 and ip from the Super-Kamiokande ratio of ratios 
in the 3rd L/E bin (see text), from the the LSND appearance signal (grey-shaded area), and from 
the Bugey disappearance measurement (hatched area left of f) = 5.5°). 

ratios at the two bins) one can fit the three mixing angles, along with an overall 
flux normalization constant. This procedure determines the mixing angles to be 
6 = 8.3°, f3 = 52.6°, and ip = 49.8°, for an overall scale factor of 0.88. The ratio 
of the measured to predicted solar neutrino flux, which is nothing else but the v e 
survival probability in Eq.(|ll|) with all sine squared terms averaged to 1/2, 

flfi- = 1 _ 2c 2 (c 2 s 2 c 2 +s 2 )) (16) 

yields Rf olar = 0.53 for this particular set of angles. However, while this solution 
indicates a significant solar neutrino deficit, consistent with the current measure- 
ments, it is not compatible with other neutrino experiments, which we discuss in 
the following sections. 
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4. The LSND Experiment 

The detector was located at 30 m from the neutrino source, with average neutrino 
energies of approximately 40 MeV and 110 MeV, for the decay-at-rest (DAR) and 
decay-in-flight (DIF) beams, respectively. We shall only consider the DAR results, 
as they are more restrictive and also more statistically significant than the DIF 
results. The — > v e oscillation probability is given by Eq.([To|). For Am 2 = 
O(10 -3 ) eV 2 the first sine squared term is negligible, while the remaining two can 
be combined by approximating AM 2 — Am 2 « AM 2 to yield 



pLSND = gin 2 {w) ^2 ^ gin 2 ( AM 2 - ) (17) 



Comparing this with the two-generation neutrino oscillation expression 

I% SND = sin 2 (2Q LSND ) sin 2 (l.27 AM 2 |) , (18) 

in terms of which the LSND event excess has been analyzed, one can effectively 
identify sm 2 (2&LSND) with sin 2 (2/3) sin 2 ip. The allowed regions obtained by this 
experiment from a preliminary analysis of the entire 1993-1998 DAR datat£l is 
shown in Fig. ||(a). The excluded regions determined from the BNL-E776LJ and 
KARMEN-2l3 experiments are also shown in the same figure. Considering these 
limits, we impose an upper bound for the AAf 2 mass scale, AAf 2 < 1.9 eV 2 , 
and a lower bound for the effective mixing, sin (2/3) sin 2 V> > 0.0014. Further- 
more, following the weak lower bound imposed by the oscillations averaging re- 
quirement, AM 2 > 0.27 eV 2 , implies an upper bound for the effective mixing, 
sin 2 (2/3) sin 2 ip < 0.086, as dictated by the boundary between the allowed region 
from LSND and the excluded region from BNL-E776. These limits determine the 
allowed region in the (/?, ip) mixing angles space illustrated by the grey-shaded area 
in Fig. |l|. The recent KARMEN-2 limit would naturally impose a more stringent 
limit of sin 2 (2/3) sin 2 tp < 0.051. However, it turns out that the restrictions imposed 
by LSND on the angles /3 and ij; are rather weak as compared to the severe limits 
imposed by the Bugey reactor experiment, which we discuss in the next section. 

5. The Reactor Experiments 

The survival probability for electron (anti)neutrinos - as typically measured by all 
reactor experiments - is given by Eq.(|ll]). The most restrictive reactor experi- 
ment for the relevant range of parameters considered in this paper is the Bugey 
experiment. EM It was carried out with three detectors at distances of 15 m, 40 m, 
and 95 m from the reactor core, with a mean neutrino energy of approximately 
3 MeV. Following the same arguments as in the LSND discussion above, the first 
sine squared term in Eq.(|ll]) is negligible, and the remaining two are combined to 
yield 

pBugey = j _ ^2 ^ ^2 ^ 2 7 AM 2 . (19) 
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Fig. 2. (a) Allowed regions in the (sin 2 (2/3) sin 2 ip, AM 2 ) space as obtained by the LSND DAR 
v e appearance experiment (grey-shaded areas). The excluded regions obtained by the BNL-E776 
and KARMEN-2 experiments are also shown (thick dashed and thick solid curves, respectively). 
The dot-dashed lines correspond to sin 2 (2/3) sin 2 ip = 0.0014 and sin 2 (2/3) sin 2 ip = 0.086. (b) 
Excluded region in the (sin 2 (2/3), AM 2 ) space as obtained by the Bugey P e disappearance reactor 
experiment. The dot-dashed line corresponds to sin 2 (2/3) = 0.037. In both (a) and (b) the dotted 
lines correspond to the AM 2 limits considered here, while the point favored by this analysis is 
indicated by the dark star. 
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A direct comparison of Eq. (|T^) with Eq. ( |l7| ) shows explicitly the difference between 
the disappearance and appearance probabilities in the three-generations neutrino 
oscillations formalism, as we have argued in the previous section. However, in 
this particular case, the expression given by Eq.(|l9|) is fully equivalent with the v e 
disappearance probability that one would write in a two-generations mixing model, 

pBugey = 1 _ sin 2 ( 2 Q Bugey ) sin 2 (l.27 AM 2 . (20) 

Therefore, one can effectively identify sin 2 (20B„ sel/ ) with sin 2 (2/3), or simply Q Bugey 
with j3. Consequently, the same limits determined for O Bugey apply to the mixing 
angle (5 in our formalism, as shown in Fig. |^(b). 

Notice that the excluded region obtained by this experiment is not shown on the 
same figure with LSND, as one typically tends to do - see for example Refs. [18,20]. 
This is because LSND, BNL-E776, and KARMEN-2 are all appearance experiments, 
whereas Bugey is a disappearance one. In the two-generations neutrino oscilla- 
tions formalism, plotting the excluded regions from a disappearance experiment on 
the same figure showing the excluded (or allowed) regions from an appearance ex- 
periment is perfectly legitimate. However, this is incorrect and misleading in the 
full thrcc-gcncrations mixing, where appearance and disappearance probabilities are 
characterized by different combinations of the mixing angles. 

For the AM 2 mass scale discussed in this Letter, the Bugey reactor experiment 
severely restricts sin 2 26 Bugey to be < 0.037, and herewith (3 — 9 Bugey < 5.5° - as 
indicated by the hatched area in Fig. |] - at the 90% confidence level. This limit is 
obviously much more restrictive than the ones obtained from the LSND discussion 
above, limiting the mixing angles f3 and ip to the two dark-shaded areas in Fig. [|. 
However, the relevance of the LSND result should not be underestimated. On one 
hand, it does not allow for to become arbitrarily small, while on the other hand - 
and more importantly - it establishes the existence of the AM 2 = 0(1) eV 2 mass 
scale, which is crucial to our approach. 

Before we go on and discuss the mixing matrix determined under these con- 
straints, we would like to first discuss a general feature that emerges directly from 
the narrow intervals to which (3 and ip are confined. Both areas of allowed values 
for f3 and if (the upper and lower dark-shaded areas in Fig. |], henceforth referred 
to as A and B, respectively) are characterized by a very low mixing angle [3 < 5.5°. 
This in turn implies that sp < 0.096 and herewith, the expression for the ratio of 
measured to predicted solar neutrino flux given by Eq.(|l6|) becomes 

R Soiar = 1-1 sin 2 (26>) + O(10- 2 ). (21) 

Consequently, a mixing angle 9 of approximately 45° will automatically imply a solar 
neutrino ratio Rf olar ~ 0.5. And this is actually exactly what happens. For values 
of f3 and ip restricted to area B in Fig. [l[ the angle 9 must satisfy 20° < 9 < 64° for 
the predicted ratio of ratios in the 8th L/E bin, as a function of 9, to agree with the 
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value measured by Super-Kamiokande at the ±lcr level, / R e = 0.456 ± 0.069. 
(For mixing angles /3 and tp in area A of Fig. [j] there is no 9 solution at the ±2tr 
level.) This implies that the solar neutrino ratio must be 0.49 < R^ olar < 0.78. 

(8) (8) 

However, imposing a stricter agreement with the central value of R^ J R\ = 0.456, 
yields values of 9 of approximately 35° and 50°, which in turn yield Rf olar ~ 0.56 
and i?f°' ar ~ 0.52, respectively. 

6. The Mixing Matrix and the LSND Mass Scale 

Fitting the ratio of ratios and the /i-like event ratios at the two L/E bins, one 
obtains the following set of mixing angles: 9 — 35.1°, (3 = 5.5°, and tp = 23.3°, 
for an overall flux normalization factor of 0.87. Herewith, the mixing matrix reads 
explicitly 

/ 0.815 0.572 0.097 \ 
U = -0.559 0.730 0.394 , (22) 
\ 0.155 -0.375 0.914 / 

and the ratio of measured to expected solar neutrino flux is predicted to be Rf olar = 
0.55. Notice that, instead of using the /i-like events, one may very well use the e-like 
events in determining the mixing angles. Despite the somewhat larger errors in the 
e-like event ratios, the result does not differ much from the previous one: 9 = 32.5°, 
(3 = 5.5°, tp = 24.6°, and thus R s e olar = 0.58. 

For the mixing angles determined here, svo 2 {2Qt i snd) = 0.0058, and thus AM 2 
is restricted to 0.5 eV 2 < AM 2 < 0.8 eV 2 , as one can easily see from the LSND 
allowed region in Fig. ||(a), below the area excluded by KARMEN-2. Thus, the 
AM 2 mass scale range favored by this analysis is well above the weak limit of 
approximately 0.27 eV 2 imposed by the oscillations averaging argument. However, 
there is yet another mass limit that has to be considered at this point, namely 
as imposed by the CDHSW disappearance experimentc 2 ] (illustrated in Fig. ||). 
This experiment was performed using a neutrino beam of approximately 2 GeV and 
two detectors at 130 m and 885 m from the mean neutrino production point. The 

survival probability is given by Eq.(^), which becomes 

pCDHSW = 1 _ 4 cos 2 p sm 2 ^ ^ _ cog 2 p gin 2 ^ ^2 ^2 ^ ^ (23) 

when neglecting the Am 2 contribution, as discussed in the previous two sections. 
Comparing this with the two-generations probability 

pCDHSW = 1 _ gin 2 {2 q cdhsw) sin 2 ^1.27 AM 2 I) , (24) 

and using the mixing angles obtained by this analysis ((3 — 5.5° and tp — 23.3°), one 
obtains sin 2 (2QcDHSw) = 0.52. For this value of the effective mixing, the values 
of AM 2 below 0.4 eV 2 are excluded at the 90% confidence level by this experiment. 
Therefore, we conclude that this mass scale must be restricted to 0.4 eV 2 < AM 2 < 



10 /. Stancu 




0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 

4 cos 2 p sin 2 \|r (1 - cos 2 p sin 2 y) 



Fig. 3. Excluded region in the (4 cos /3sin 2 ^(l — cos /3sin 2 ip), AM 2 ) space obtained by the 
CDHSW Ufj, disappearance experiment, along with the point favored by this analysis (dark star). 

0.5 eV 2 , and we set it to AM 2 ~ 0.45 eV 2 for the discussion below. With these 
parameters, the point favored by this analysis lies at the edge of the LSND allowed 
region and also at the edge of the Bugey and CDHSW excluded regions, as shown 
in Figs. 0(a), Bb), and 0, respectively. 

With the mixing matrix fully determined and with Am 2 = O10~ 3 eV 2 and 
AM 2 ~ 0.45 eV 2 , we can also investigate the atmospheric neutrino ratio of ratios 
predicted by this analysis as compared to the Super-Kamiokande data. This is 
illustrated in Fig. [|(a) for Am 2 = 10~ 3 eV 2 , as we shall make no attempt in this 
Letter to fit this mass scale. Our predicted /x-like and e-like event ratios are shown 
in Fig. ||(b), together with the Super-Kamiokande data rescaled by an overall factor 
of 0.87. In both cases the agreement between our prediction and the experimental 
data is reasonably good. No smearing of L/E has been included in our simulations. 

7. Conclusions 

We have shown that the L/E dependence of the Super-Kamiokande atmospheric 
neutrino data fully determines the three-generations neutrino mixing matrix, as 
given by Eq.(|2^). The only underlying assumption is that the mass scales involved 
are those determined by the atmospheric neutrino data, Am 2 = O(10 -3 ) eV 2 , and 
by the LSND experiment, AM 2 = 0(1) eV 2 . Consequently, the ratio of measured 
to predicted solar neutrino flux is determined to be Rf olar — 0.55, independent 
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Fig. 4. (a) The predicted ratio of ratios (i? M /i? e ) versus log(L/ E) - dashed histogram - compared 
to the Supcr-Kamiokande data (empty crosses with error bars), (b) The predicted /i-likc (dotted) 
and e-like (dashed) event ratios (i? M and R e ) versus log(L/E), compared to the Super-Kamiokandc 
data (full and empty circles with error bars, respectively) scaled by a factor of 0.87. Both (a) and 
(b) have been obtained with Am 2 = 10~ 3 eV 2 and AM 2 = 0.45 eV 2 . 
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of energy. While no attempt has been made here to determine a more precise 
value for Am 2 , we believe that AM 2 ~ 0.45 eV 2 . This, together with the effec- 
tive sin 2 (29iSArc) = 0.0058 determined here, lies just outside the sensitivity of 
the KARMEN-2 experiment, to be reached at the end of 2001. However, the Mini- 
BooNE experiment at FermiLab should see a significant number of — > v e events at 
this mixing: approximately 345 (525) per year for AM 2 = 0.4 (0.5) eV 2 . The long- 
baseline — * v T appearance experiments (K2K, MINOS), which are probing the 
Am 2 mass sale, should also see significant excess events for this particular neutrino 
mixing pattern. Higher statistics in the Super-Kamiokande atmospheric neutrino 
data and a more detailed analysis - which may include deviations of r = N^/N e 
from the naive value of r = 2 due to geo- magnetic effects and L/E variations - 
should further restrict the mixing angles in this neutrino oscillations model, which 
involves the Am 2 = 0(1O~ 3 ) eV 2 and AM 2 = 0(1) eV 2 mass scales. However, 
should the LSND result prove to be wrong, it might be just as simple as: atmo- 
spheric deficit = Ufj, — ► u T , and solar deficit = v e — » v x . On the other hand, should 
the LSND result stand and should the solar neutrino deficit show a strong energy 
dependence, more complex neutrino models have to be invoked. Either way, if neu- 
trino oscillations are firmly established, this will open exciting avenues beyond the 
physics of the Standard Model. 
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